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The superconducting order parameter and electron-phonon interaction (EPI) in a MgB2 c-axis 
oriented thin film are investigated by point contacts (PCs). Two superconducting gaps Asmaii ~ 
2.6 meV and Ai nrge ~ 7.4 meV have been established, which fairly well correspond to the theoretical 
gap ratio f :3. The reconstructed PC EPI function along c-axis has above 30 meV similarities with 
the phonon DOS, while a pronounced maximum around 20 meV is resolved pointing out that the 
low frequency phonon mode can be of importance for pairing. We have estimated the Fermi velocity 
in MgE>2 about 5xl0 7 cm/s. 
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INTRODUCTION 

The compound MgB2 with graphite-like planes of 
boron atoms attracts much attention as superconduc- 
tor with at the present the highest T c ~40K |l| for 
a binary systems. The observation of boron isotope 
effect [2, 3 1 and examination of electron-phonon cou- 
pling [I, 5, ||, 0, H [| [HJ in MgB 2 are in accor- 
dance with the expectations for conventional BCS su- 
perconductivity mediated by electron-phonon interaction 
(EPf). Investigations of the order parameter by tunncl- 

and point-contact tech- 
22J, |23j |2j confirm that MgB 2 is 
most likely an s-wave superconductor. In all cases the 
spectra show unambiguous features of an energy gap A 
in the density of states (DOS), albeit the results are con- 
troversial on the gap width. Values of A ranging from 
1.5 to 8meV have been reported, pointing out the possi- 
bility of an anisotropic or distributed (nonhomogeneous) 
energy gap or even multiple gaps. The latter scenario 
has been recently recalled by Liu et al. || for MgB 2 . 
The Fermi surface consisting from nearly cylindrical hole 
sheets arising from quasi-2D boron bands and three di- 
mensional tubular network Q was considered. The dif- 
ferent character of the sheets raises the possibility that 
each has a distinct gap. The existence of two different 
energy gaps with the ratio approximately 1:3 being re- 
spectively smaller (for the 3D gap) and larger (for the 
2D gap) than the standard weak coupling BCS value 
A=1.76k B T c ~6meV was predicted in @. 



Determination of the Eliashbcrg function a 2 F for su- 
perconducting systems provides a consistency check of 
the phonon-mediated pairing mechanism. In j^, |^, [t], 
§, g, the phonon dispersion and DOS, the EPI 
function a 2 F and electron-phonon coupling constant A, 
T c and isotope effect for MgB 2 have been calculated 
using different methods and approximation. Estima- 
tion of the electron-phonon-coupling strength in MgB 2 
@ 1 i E E i y ields A - L Established by different 
authors, the phonon DOS and the EPI function have 
both similarities and differences. The common feature 
for presented in [pi pL Oil, 25 phonon DOS is the maxi- 
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mum energy of about 100 meV, absence of visible phonon 
peaks below 30 meV and rich structure with a number of 
peaks in between 30-100 meV. The main feature of the 
calculated EPI function is a dominant maximum between 
60 and 75 meV arising from in-plain boron E 2g phonon 
mode. The precise energy position and shape of the max- 
imum depend on the method of calculation || |(], §j] . The 
transport EPI function, as shown in Q , mimics the EPI 
one, although, as mentioned in ||, interband anisotropy 
reduces the transport coupling constant from 1 to 0.5. 
From both the DC resistivity and optical conductivity 
measurements |26 , even a smaller value of X tr = 0.13 is 
derived. 

However, some neutron inelastic scattering measure- 
ments with higher resolution J27], |2^] revealed a clear 
maxima in generalized phonon DOS below 30 meV, 
namely at 16 and 24 meV in |27| or at 17.5 meV in p8[ . 
Therefore it is an open question whether the only high 
frequency phonons are responsible for thermodynamic, 
transport and especially for fascinating SC properties of 
MgB 2 . 



It has been also suggested in 1 29 that the " multi-gap" 
structure observed in the tunneling spectra fll5| can be 
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explained by considering of a low-frequency phonon mode 
at 17.5 meV, which is revealed in the inelastic neutron 
scattering experiment | p8[ . Soft bosonic mode was also 
exploiting in ]3C| to describe an upper critical field be- 
havior of MgB2 within a multi-band Eliashberg model. 

Measuring of a nonlinear conductivity of point contacts 
(PC) between two metals allows in a direct way recover- 
ing the EPI function a 2 F pl[, a s well as investigating the 
superconducting (SC) gap ]32[. The goal of this report 
is simultaneous study of peculiarities of EPI along with 
SC order parameter of MgB2 by PC method in order to 
clarify the mentioned issues, namely, the multigap struc- 
ture of the order parameter and the peculiarities of EPI 
function, trying to find their correlation. 



EXPERIMENTAL AND RESULTS 

We have used the high-quality c-axis oriented MgPj2 
thin films grown by a pulsed laser deposition technique 
p3[ . The films were grown on (1 10 2) AI2O3 substrates. 
The typical film thickness was 0.4 Ltm. The resistivity of 
the film exhibits a sharp transition at 39 K with a width 
of ~ 0.2 K from 90% to 10% of the normal state resistivity 
|Q . The residual resistivity po at 40 K is ~ 6 liQ, cm J3EJ 
and RRR=2.3. 

Different contacts were established in situ at liquid 
4 He temperatures by touching as prepared surface of the 
MgB2 film by the sharpened edge of an Ag counterelec- 
trode, which were cleaned by chemical polishing. By this 
means the current flow through PCs is preferably along 
the c-axis. The experimental cell with the samples holder 
was immersed directly in liquid He to ensure good ther- 
mal coupling. Both the differential resistance dV/ dl and 
d 2 V / dl 2 (V) vs V were registered using a standard lock- 
in technique. The resistance Rn at y > A of investi- 
gated contacts ranges from 10 to 1000 fl at 4.2 K. 



Superconducting gap 

According to the Blonder-Tinkham-Klapwijk theory 
p2| of conductivity of N-c-S metallic junctions (here N is 
normal metal, c is constriction and S is superconductor) 
a maximum at zero-bias voltage and a double-minimum 
structure around V ~ ±A/e in the dV/dl curves man- 
ifest the Andreev reflection process with a finite barrier 
strength parameter Z. Thus the positions of the minima 
roughly reflect the SC gap value. This follows from the 
equations for the current-voltage characteristics 
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FIG. 1: Reduced differential resistance R^dV/dl vs. V mea- 
sured for three MgB2-Ag contacts with one gap minima (curve 
1, Rn=350 Q) and two gaps (curve 2,3, -Rjv=20 and ~25f2) 
structures. Symbol traces and bottom solid line are the ex- 
perimental curves at T = 4.2 K. Thin dashed lines are theo- 
retical dependence calculated using Eq.(l) with one gap 2.7 
meV (curve 1, Z—0.9, r=1.2 meV) and two gaps: 7.1 and 2.7 
meV (curve 2, Z=0.9, T=0.26 meV), 7.4 and 2.6 meV (curve 
3, Z=0.64 and 0.38, T ~ 0). The gap weight factor ratio 
W small /W Large is 10 and 1.5 for curves 2 and 3, respectively. 
The y-axis one division corresponds to 0.05 for curve 1 and 
0.1 for curve 2. 
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where f(e) is the Fermi distribution function. The broad- 
ening of the quasiparticle DOS in the superconductor can 
be taken into account if e will be replaced by e — iT in Eq. 
(1). We used Eq.(l) to fit the measured dV/dl curves of 
PCs and extract SC gap. 

The dV/dl curves taken at 4.2 K <C T c show usually 
structure, which contains symmetric with respect to zero 
bias minima at about ±2-3 mV with zero-bias maximum 
at V=0 (see Fig. 1, curve 1), related to the SC gap A. 
Often the dV/dl vs. V curves show clear two gap struc- 
ture with shallow features corresponding to a larger gap, 
as shown by curve 2 in Fig. 1. Among the many dV/dl 
curves (about 100), there was a curve 3 (Fig. 1) exhibit- 
ing almost equal features corresponding to the both gaps 
with minima at about ±7 and ±2.3 mV. These values 
are close to the maxima at 7 and 1.5-2 mV of the minima 
distribution in dV/dl given in pjJ. Fitting of curve 3 
by Eq. (1) yields the gap value 7.4 and 2.6 meV, which is 
consistent with PC measurements 7 and 2.8 me V in |2l[] 
or 6.2±0.7 and 2.3±0.3meV in j24|] and with single small 
gap 2.6 meV p3|. In our case, the ratio 2.85 of the larger 
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gap to the lower one is close to the theoretical value 3:1 



Phonons features 

The voltage V applied to the ballistic contact defines 
the energy scale eV for scattering processes. EPI re- 
sults in backflow scattering of energized electrons so that 
some of them are reflected by creating of phonons and do 
not contribute to the current. This leads to a nonlinear 
I — V characteristic and, as shown by the theory of Ku- 
lik, Omelyanchouk and Shekhter |36| , second derivative 
d 2 V/dI 2 (V) of the I — V curve is proportional at low 
temperatures to ap C F(uj): 
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where ape, roughly speaking, describes the strength 
of the electron interaction with one or another phonon 
branch and takes into account the kinematic restriction 
of electron scattering processes in the point contact by 
factor 1/2(1 — #/ tan#) (for transport and Eliashberg EPI 
functions the corresponding factors are (1 — cos£?) and 1, 
respectively), where 9 is the angle between initial and 
final momenta of scattered electrons. It should be men- 
tioned that in PC the large angle 9 — > tt scattering (back- 
scattering) processes of electrons are strongly underlined. 

The PC resistance is determined by a sum of ballis- 
tic Sharvin and diffusive Maxwell term according to the 
simple formula derived by Wexler [B8| : 
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which is commonly used to estimate the PC diameter d. 
Here pi = p^/ne 2 , where pp is the Fermi momentum, n 
is the density of charge carriers. The latter for MgB2 is 
estimated in n ~ 6.7 x 10 22 |37j], which results in pi ~ 
2 x 10~ 12 f2-cm 2 using vp ~ 5 x 10 7 cm/s Q. Hence, 
the upper limit for elastic mean free path for our film is 
about 3 nm. In this case, using Eq. (3), we can find that 
condition d < I is fulfilled for PC with R > 40 £1 or fol- 
lower resistance supposing multiple contacts in parallel. 

The measurements of d 2 V/dI 2 (V) dependencies to 
recover spectral EPI function reveal a wide variety of 
curves. Nevertheless, we were able to select similar 
d 2 y/d/ 2 (y) characteristics for PC with different resis- 
tance (see Fig. 2). The common features are the repro- 
ducibility of position of the main maxima placed at about 
40, 60 and 80-90 meV and lack of spectral features above 
100 meV. There is a correlation in the peak position in 
d 2 V/dI 2 (V) and calculated in § EPI function (Fig. 2). 

The mentioned features are also seen in the spectrum of 
another PC in Fig. 3. Note, that this curve is measured in 
a magnetic field of 4T, which have to suppress additional 
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FIG. 2: Second derivative d 2 V /dI 2 (V) averaged for opposite 
polarities for 4 MgB2-Ag contacts with different resistance 
(R = 43, 36, 111 and 45 fi from top curve to bottom one) at 
T=4.2 K with clear maxima in the range of phonons above 30 
meV. Bottom curve is calculated in || smoothed with exper- 
imental resolution (the segment above) EPI function. Signal 
increase in experimental curves below 30 meV is due to the 
SC gap structure. 



features on d 2 V / dl 2 (V) , which could arise from SC weak 
links or degraded SC regions. It is also worth to be noted, 
that dV/dI(V) of this contact increases above 30 mV, in- 
dicating a direct metallic contact. Absence of a barrier 
layer at the interface is confirmed by proximity induced 
superconductivity resulting in a sharp dip at V—0 on 
dV/dI(V) (see Fig. 3, inset). Therefore, it is also reason- 
ably to suggest for this contact that the barrier parameter 
Z, obtained from the fit (see Fig. 3, inset), is caused only 
by the mismatch of the Fermi velocities between the two 
electrodes. According to 1 32 Z = (1 — r)/2r 1 / 2 , where 
r is the ratio of the Fermi velocities. Using the fitting 



parameter Z=0.5 for this contact we find the Fermi ve- 
locity of MgB 2 about 5xl0 7 cm/s taking 1.4x10 s cm/s 
as the Fermi velocity of Ag. 

After subtracting a linear background, as shown in 
Fig. 3, the presented in Fig. 4 PC EPI spectrum (curves 
1,2), has above 30 meV similarities with the phonon DOS 
measured by neutron scattering ]9], [n|. Below 30 meV 
all the spectra (see Fig. 2, 3) exhibit a steep increase con- 
nected with superconductivity (Andreev reflection). We 
have tentatively subtracted the nonlinear background in 
this region, as shown in inset of Fig. 4, and found an in- 
tensive peak at 20 meV and around 30 meV (Fig. 4). In 
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FIG. 3: First dV/dJ(V) (dashed) and second d 2 V/dI 2 (V) 
(solid) derivative of I — V characteristics for MgB2-Ag con- 
tact at 4.2 K and magnetic field 4 T. The d 2 V/dI 2 (V) curve 
is truncated near zero bias due to large signal caused by the 
SC gap. Straight line shows an approximate background cor- 
responding to the parabolic shape of the first derivative. In- 
set demonstrates fitting of dV/d/(V) (symbols connected by 
dashed line) by Eq.(I) (thin solid line). Parameters of the 
fit are: A=3.8 meV, Z=0.5, F=3 meV. Dip in dV/dI(V) at 
V=0 is tentatively due to the proximity effect. 



principle, these peaks (as the other lower-frequency peaks 
not revealed yet) might not related directly to the EPI 
function, and could be due, e.g., to the suppression of the 
SC order parameter. This suppression can occur near the 
characteristic phonon energies with small group veloc- 
ity, corresponding to peaks in phonon DOS, as observed 
for PC spectra of ordinary superconductor tantalum p9[ . 
Indeed, corresponding peaks may be found in the men- 
tioned energy region at 16-17, 24, 31 meV [^7], or 
hillock at 20 meV (^] in the phonon DOS determined by 
neutron scattering. Interesting, that PC EPI spectrum 
of pure Mg has maxima at 17 and 28 meV (see Fig. 4a) 

From the comparison of our data with the phonon 
DOS follows that all phonons above 30 meV contribute 
in average with approximately equal weight to the PC 
EPI function. We do not see the prevailing of one mode 
around 60-70 meV, as follows from calculated Eliashberg 
or transport EPI function in Ref.|^, || £|. Probably, the 
reason is that our measurements are mainly along the 
c-axis. 

Concerning the estimation of an absolute value of EPI 
function or EPI constant A from PC data, it is prema- 
turely because the nonlinearity of PC I—V curves in the 
region of characteristic phonon frequencies is too low. As 
is seen from Fig. 3, dV/dl increases above 30 meV only 
on a few percent contrary to 10-50% increase in the zero 
bias region due to the Andreev reflection at zero mag- 
netic field. The short elastic electron mean free path 
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FIG. 4: a) PC EPI function (see Eq.(2)) of MgB 2 (curves 1 
and 2) reconstructed from d 2 V /dl 2 (V) in Fig. 3 averaging for 
plus and minus polarity after subtracting linear (as is shown 
in Fig. 3) and hand-made (see inset) b ackg round. Solid curve 
(3) shows PC spectrum of pure Mg |4C|| . Inset: proposed 
nonlinear background (symbols) below 30 meV for the same 
contact recorded with enlarged scale, b) Phonon DOS accord- 
ing neutron measurements || (solid line) , (symbols) and 
j27[ (symbols connected by dashed line). In the latter case 
only low energy (<40meV) part of phonon DOS is shown. 



in the constriction (diffusive regime) is likely the main 
reason of the small nonlinearity. Sure, to receive both 
quantitative results and final shape of EPI function the 
investigations of more perfect samples or even single crys- 
tals are very desirable. 



CONCLUSION 

Our results on MgB2 c-axis oriented thin film investi- 
gated by point contacts unequivocally indicate the pres- 
ence of two SC gap with the theoretically predicted ratio 
1:3 J8j. Above SC gap the EPI features dominate in the 
PC spectra. Thus the electron-phonon coupling in MgB 2 
must therefore be included in any microscopic theory of 
superconductivity. The PC EPI function has similarities 
above 30 meV with the phonon DOS, while maxima at 20 
and 30 meV are also resolved pointing that these low fre- 
quency phonon modes should be considered by describing 
of SC state of MgB2. Using the barrierless junction, the 
Fermi velocity 5x 10 7 cm/s in MgB2 have been estimated. 
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